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Time-dependent nonequilibrium properties of a strongly correlated electron system driven by 
large electric fields is obtained by means of solving the time-dependent Schrodinger equation for the 
many-body wave function numerically in one dimension. While the insulator-to-metal transition 
depends on the electric field and the interaction, the metallization is found to be described in terms 
of an universal Landau-Zener quantum tunneling among the many-body levels. These processes 
induces current oscillation for small systems, while give rise to finite resistivity through dissipation 
for larger systems/on longer time scales. 
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Introduction Nonlinear responses and time-dependent 
phenomena are an open frontier in the physics of strongly 
correlated systems. Namely, while the physics of elec- 
tron correlation has been intensively studied in the last 
decades in the context of the high-Tc cuprates and re- 
lated systems, the efforts were concentrated on equilib- 
rium or linear-response properties. Properties beyond 
the linear response or time-dependent properties have 
been much discussed] 1], but remains a challenging prob- 
lem, which is, theoretically, due to the difficulty in deal- 
ing simultaneously with the many-body effect (correla- 
tion) and finite current (non-equilibrium), which are both 
non-pert urbative physics. So we definitely require a clear 
picture for this fundamental problem, especially for low- 
dimensional materials. 

In this Letter we consider the Hubbard model un- 
der constant driving forces in one dimension. The 
ground state of the half-filled Hubbard system is a 
Mott insulator for arbitrary strengths of the electron- 
electron repulsion /7 > 0, while the state is metallic when 
the band filling is shifted away from the half-filling by 
doping carriers. So the question we pursuit here is: what 
will happen if we destroy the Mott insulator by applying 
strong electric field, instead of doping? 

There are existing theoretical approaches that employ 
the Bethe-ansatz method, and the closing of the Mott 
gap has been discussed j^, H- In these studies, how- 
ever, electric fields are not actually applied, but the left- 
going and right-going hopping terms are made different 
instead. It is rather difficult to relate this artificial, and 
non-Hermitian, model with a system in an electric field. 
That is why we have here opted for actually applying 
electric fields for the first time to keep track of the time 
evolution of the many-body wave functions and levels. 
For that we make use of a numerical integration of the 
time-dependent Schrodinger's equation, where we have 
also looked at excited states, which should be essential 
in examining non-equilibrium phenomena. 

We have found that the electric field, if strong enough, 
breaks the Mott-insulator phase. While the critical field 
strength required for the breakdown of the Mott insula- 
tor depends sensitively on the magnitude of the electron- 
electron interaction, we propose here that the mecha- 



nism for the metallization can be viewed as the non- 
adiabatic tunneling between the many-body levels. We 
have verified this by confirming numerically that a uni- 
versal Landau-Zener quantum tunneling governs the non- 
linear conduction (e.g., the I — V characteristics repre- 
senting the insulator-to- metal transition), where the fit is 
surprisingly good although the Landau-Zener formalism 
is originally intended for one-body problems while the 
problem at hand, being many-body, involves the Hilbert 
space with huge dimensions. 

Just after the metallization, we have a self-induced 
current oscillation. While this should be realistic for 
mesoscopic systems, for large systems or on a longer 
time scale, a novel, "Ohmic" conduction is found to re- 
sult. This occurs despite the absence of disorder and 
the heatbath degrees of freedom, but a series of non- 
adiabatic tunneling among many-body states are respon- 
sible for the dissipation effect. Indeed, the breakdown of 
a one-dimensional (ID) _5] as well as two-dimensional 
(2D) 6] Mott insulator has been experimentally studied, 
and, among other interesting phenomena including spon- 
taneous density-pattern formations,!^ 0, a seemingly 
Ohmic conduction was found for a rather wide range of 
external electric field, until the Ohmic conduction is even- 
tually broken. 

Formulation If one wants to apply a finite electric field 
to a system, one problem is how to treat the electrodes. 
To get rid of ambiguities arising from this, we have opted 
here for a periodic system (a ring in ID), where the 
electric field F is applied via a time-dependent AB flux 
^(t) = eLFt piercing the ring (inset of Fig^a); L: sam- 
ple length). This will lead to a circular electromotive 
force due to Faraday's law. The flux makes the hop- 
ping integral in a tight-binding model complex, where 
the Hamiltonian is 
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Here W is the bandwidth, U the electron-electron repul- 
sion, N = L/a the total number of sites with a being 
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the lattice constant. Hereafter we take the unit in which 
e = a = h = 1. The electric field in this formalism 
amounts to a time-dependent deformation of the Hamil- 
tonian. 

We have then to solve the time-dependent Schrodinger 
equation, = which governs time 

evolution of the quantum system at absolute zero tem- 
perature, starting from the ground state of the Hamil- 
tonian at t = 0, = 0)) = |^o), which is obtained 
with the Lanczos method. The time-integration of the 
state vector, which, being many-body, has a huge di- 
mension, requires a reliable algorithm. So we adopt 
here Cranck-Nicholson's method that gurantees the 
unitary time evolution, where the time evolution is put 
in a form, \^{t ^ At)) = e-'I'^^'^^'^ dt \^{t)) - 
[1 - iAt/2H{t + At/2)]/[l + iAt/2H{t + At/2)] 
which is unitary by definition. Here the time step is taken 
to be small enough {dt = 1.0 x 10~^ with the time in units 
of Ah/W hereafter) to ensure convergence for < 10 site 
systems, for which the dimension of the Hamiltonian is 
~ 10^. We have concentrated on the total = ^ sub- 
space with = Ni = N/2. 

Result We first plot in Fig^a) the result for the ex- 
pectation value of the current density averaged over the 
sites, 
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The behavior of J{t) for various values of U with a fixed 
value of the electric field F is seen to fall upon three 
regimes when U is varied: A perfect metallic behavior 
(J(t) oc t) when the electrons are free {U /W = 0), while 
when the interaction is strong enough (U /W ^ 1) the 
current has a zero expectation value. For an intermediate 
regime of U /W we have finite J's with some oscillations 
in the current for finite systems. By contrast, a non- 
half-filled system has the time evolution distinct from 
the ground- state behavior. The difference has its root in 
the spectral property as will be discussed later. 

Now, FigH^b) plots the time evolution of the current 
when the electric field F is varied with a fixed value of 
again for the half-filled case. The result may be 
summarized as follows: 

(i) Small F regime (Mott insulator) 

A drastic difference between the half-filled and 
doped systems appears for small F. When half- 
filled, J{t) in the limit of F ^ smoothly ap- 
proaches the {J{t)) = behavior of the ground 
state (Mott insulator). Here {Jit)) is the time- 
averaged current. On top of the {J{t)) = an os- 
cillatory behavior with the period of $o(= = 1: 
fiux quantum) is seen, which is nothing but the AB 
oscillation (a saw-tooth, due to a symmetry about 
^ = ^o/2). 

(ii) Moderate F regime (metal) 

In this regime, the current in the half-filled case 
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FIG. 1: (a) Time evolution of the current, J(t), for a half- 
filled, 10-site Hubbard model for various strengths of the cor- 
relation, < IJ/W < 5, for a fixed electric field F = 1/lOL. 
Time is measured in units of rw = 4:h/W, LF in VK/(4/i), 
and J{t) in 1/rw- The range of the time in this panel corre- 
sponds to a range of the AB-flux < $/$o < 1- The inset 
schematically depicts the sample geometry, where an AB flux, 
$(t) = LFt, increasing linearly with time induces an elec- 
tric force through Faraday's law. (b) A longer-time behavior 
of the current when F is varied with a fixed U/W — 0.25, 
again for the half-filled case. Here the horizontal axis is 
LFt = $(t)/$o and over < $/$o < N{^ 10 here). The 
inset schematically shows three kinds of behavior (MI: Mott 
insulator, M: metal, PM: perfect metal), (c) A plot similar to 
(b) for a non-half-fllled case (iV^ = AT^ = 3 < N/2 = 5). 



shows an oscillatory behavior (see typically the 
LF = 0.008 data in FiglUb)). 

(iii) Large F regime (perfect metal) 

When the electric field F becomes large enough, 
the system becomes a metal, in which {J{t)) oc t 
for A<l> < <l>oA^/4. A further oscillation in J(t), 
with a long period (A^ = ^qN), is seen, which we 
will discuss later. 

The F-dependence of (J) are displayed in Fig|21as the 
FV characteristics for various values of U. Here the time 
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average (J) = J^\j{t))dt/ At {LFAt = N^o/^) of the 
current density J{t) is taken over one forth of the ex- 
tended AB period (0 < <l> < A/'<l>o/4), since we are in- 
terested in the rise in the current which should represent 
the behavior in the thermodynamic hmit. We can see 
that (J) becomes nonzero rather abruptly at the metal- 
lization, where the threshold electric field increases with 
U /W, thereby the F- dependence becoming weaker. Just 
after the metallization some oscillation (in the F- depen- 
dence this time) is seen for finite systems. 
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FIG. 2: FV characteristics for various values of U/W for the 
half- filled Hubbard model with N — 6. 




FIG. 3: (a) Landau- Zener process around a level anticrossing, 
where AE denotes the gap, p the transition probability, and 
SE the difference between the "unperturbed", crossing energy 
levels (dashed lines), (b) [/-dependence of the many-body 
gap marked with a double circle in (c). (c) The low- lying 
levels versus $ for the Hamiltonian eq.(0 in the half- filled 
case (N = 10 with = = 5) for U/W = 0.125. Level 
repulsions due to the interaction U are encircled, (d) Similar 
plot for a doped case {N = 10 with = Ni = 3). 



Non-adiabatic tunneling In order to understand the 
physics underlying these time evolutions and their F de- 
pendence, let us here evoke the notion of non-adiabatic 
tunneling, originally conceived for one-bod y pr oblems by 
Landau, Zener, and Stiickelberg fLZS)[^. llOl ITU . The 
LZS theory considers the quantum tunneling beyond the 
adiabatic approximation. Namely, when a parameter 
defining the Hamiltonian is varied infinitely slowly (adi- 
abatically), we can just plot a set of energy levels against 
the parameter, which in general contain level anticross- 
ings, since the two levels should repel with each other 
unless they are allowed to cross due to e.g. symmetry 
reasons. An initial state that start from one of the lines 
should evolve with time sticking to that line (adiabatic 
theorem, corresponding to p = in Fig|3fa)). When the 
parameter is varied with a finite velocity, the state can 
make a transition across the level anticrossing with a fi- 
nite probability p{y^ in FiglHIJa)), where the transition 
is caused by a quantum mechanical tunneling across the 
gap AE. The transition probability p depends in the LZS 
theory on the speed {6E) the two energy level approach 



as 
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Here 5E is the difference between the "unperturbed", 
crossing energy levels (dashed lines in Fig|3Ja)) and 
5E = dSE/dt = {d6E/d^)^ with LF = d^/dt. We 
can immediately see that the process is non-perturbative 
in F, since p is singular in F. 

Although the original LZS theory is devised for one- 
body systems, there is no reason why we cannot apply it 



to many-body systems, as demonstrated for a spin sys- 
tem by Miyashita et al[l3,[l3- So here we apply the con- 
cept to the Hubbard model, which is, to our knowledge, 
the first time the LZS theory is applied to interacting 
electron systems. In order to check the validity of the 
LZS picture, we have first calculated the transition prob- 
ability p. The level anticrossing we focus on is the first 
one encountered by the ground state at ^ = /2 in the 
level flow (marked with a double circle in Fig|3fc)). In 
FigSIa) we plot |(^o|^(^))P, the weight of the ground 
state around the level anticrossing, for various values of U 
and F. From this we have obtained the transiti on p roba- 
bility p as the asymptotic value of | (^o|^(^)) Fig- 
ure ^b) plots -Log p as a function of the LZS parameter, 
{AE)'^/{SE)^ with AE now defined as the interaction- 
originated one. We can see a remarkably accurate linear 
dependence on the LZS parameter, which clearly indi- 
cates that the LZS theory is applicable to the many- 
body system we have at hand. So the field F and the 
interaction U enter into the problem only as a combina- 
tion {AE{U)f/{^LF) (where AE{U) is an increasing 
function of /7 ; see FiglH^b)). 

For large enough electric field F the nonadiabatic tun- 
neling is so effective that the state goes straight across 
each crossing with probability close to unity, so we end 
up with a long period (A<l> = ^qN)., an analogue of the 
"extended AB period" discussed for electron systems by 
Kusakabe and two of the present authors 15, 16], a no- 
tion originally proposed for a spin (Heisenberg) system by 
Sutherland [l3- In the thermodynamic limit [N oo) 
the period becomes oo. 

Encouraged by this, we have then re-plotted the time- 
averaged current, (J), against the inverse of the LZS pa- 
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in the half-filled Hubbard model, which is a charge gap 
characteristic to the half-filled Hubbard model and van- 
ishes when doped (Figin^c)). The total momenta of 
the anticrossing states (the ground state and the ex- 
cited state wi th one pair of complex charge rapiditie s [ T^) 
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FIG. 4: (a) Time-evolution of the weight of the ground state, 
|($o|^(t))|^, calculated for the half-filled Hubbard model 
{N — 10, A^i = Ni — 5) for various values of F with 
U/W = 0.025. Inset shows the solutions [T^ of the LZS equa- 
tion with its asymptotic values indicated, (b) Transition prob- 
ability p (decrease in |($o|^(t))|^) plotted against the LZS 
parameter ^^^^ {6E = ^^LF) for various values of F and 
U/W. 

rameter ^/(AE)^ x LF, in FigEl Dramatically, ah 
of the curves, which appeared quite different for differ- 
ent values of U /W in the raw I-V characteristics (FigjSJ, 
fall on a single, universal curve within a reasonable er- 
ror when plotted against the LZS parameter. Specifi- 
cally, a threshold between the insulating behavior and 
the dissipative metallic one is clearly seen at around 
F~0.5x(A^)V(^^L). 

After the metallization the current is seen to behave 
roughly linearly with F. This seemingly "Ohmic" be- 
havior is far from trivial but rather surprising, since 
we have a many-body but clean system. This implies 
that (i) after many level crossings the system reaches a 
steady state, and (ii) the many-body gaps at these level 
crossings have similar magnitudes as the first one {AE 
above). We can in fact recall that the non-adiabatic tun- 
neling is a quantum version of dissipation, where differ- 
ent quantum states become mixed after level anticross- 
ing. To get an insight into the quantum dissipation let 
us examine the nature of the many-body gap AE{U) 



FIG. 5: Time- averaged current plotted against the inverse 
LZS parameter (dSE /d^) /{AEf x LF for various values of 
U/W and F in the half-filled Hubbard model with N = 6 
(open symbols) and N = 8 (solid ones). Inset is a blowup of 
the same data for small F on a log-scale. 



differ by 2kF where kp is the Fermi wavenumber (for 
small U ; replace kp with the quasimomentum for gen- 
eral U). At half-filling, these Umklapp processes take 
place and the momentum of the many-body state is dis- 
sipated. In other words, the role of heatbath degrees of 
freedom is played by the many-body system itself. 
Let us add that the current oscillation (typically seen for 
LF = 0.008 in Fig^b)) is due to the kicks from Umklapp 
processes and should be observed in small systems with 
strong electron correlation. The threshold electric force, 
F ^ 0.5{AE)'^ /[{d6E/d^)L], has a similar order of mag- 
nitude as the critical field observed in the experiment 
The threshold electric force, F - {AEf /[{dSE/d^)L], 
should approach an asymptotic value in the thermody- 
namic limit, since {dSE /d^)\^=^Q/2 «l>=0 ^ 
Drude weight) ^ 1/L should cancel the L in the denom- 
inator. Exactly how the thermodynamic limit is reached 
and how the seemingly Ohmic conduction occurs is an 
important future problem. 

We thank Koichi Kusakabe for discussions, and TO 
acknowledges Naoki Watanabe for a help on the time- 
dependent numerical simulation. 
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